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Non-hexagonal-ring defects and structures induced by strain in graphene and in 

functionalized graphene 

Joice da Silva-Araiijo, A. J. M. Nascimento, Helio Chacham, and R. W. Nuneqj 

Departamento de Fisica, Universidade Federal de Minas Gerats, 
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We perform ah initio calculations for the strain-induced formation of non-hexagonal-ring defects 
in graphene, graphane (planar CH), and graphenol (planar COH). We find that the simplest of 
such topological defects, the Stone- Wales defect, acts as a seed for strain-induced dissociation and 
multiplication of topological defects. Through the application of inhomogeneous deformations to 
graphene, graphane and graphenol with initially small concentrations of pentagonal and heptagonal 
rings, we obtain several novel stable structures that possess, at the same time, large concentrations 
of non-hexagonal rings (from fourfold to elevenfold) and small formation energies. 

PACS numbers: 



I. INTRODUCTION 

The intense interest in the Physics of graphene, over 
the last few years, derives in large part from the ex- 
ceptional electronic properties of the pristine mate- 
rialii"— While pristine samples obtained by exfoliation 
of graphite have allowed exploration of a rich variety of 
phenomena connected with the Dirac-fermion nature of 
the electronic excitations in graphene, current large-scale 
synthesis protocols have reached the point where poly- 
crystalline or highly-inhomogeneous samples are being 
produced. ^^- In very recent experiments, graphene layers 
resulting from reduction of graphene oxide show conduc- 
tivities that are 1 or 2 orders of magnitude smaller than 
those of pristine graphene ^Sr-i^ possibly due to residual 
dopants and intrinsic defects. Better conductivities can 
be obtained when partially oxidized graphene is used as 
a precursor ji^ 

From the perspective of applications, widespread use 
of graphene will depend on the control and understand- 
ing of the properties of such defective graphene products. 
The atomic structure of defects in reduced graphene ox- 
ide (RGO) has been recently investigated through high 
resolution transmission electron microscopy (TEM)^ii^ 
and molecular dynamics simulations. S The experimental 
results show significant areas with holes, grain bound- 
aries, isolated pentagon-heptagon pairs and clusters of 
carbon pentagons, hexagons, heptagons and a few oc- 
tagons. Similar defect morphologies can be generated in 
exfoliated graphene by electron-beam irradiation,® and 
by irradiation under extreme temperatures J The occur- 
rence of areas with high density of extended topological 
defects (ETD) observed in the samples of Ref. 5 was at- 
tributed to the process of oxidation and further reduction 
of the graphene oxide, since such large defective areas do 
not appear in mechanically exfoliated graphene from the 
same graphite source. Moreover, in these samples resid- 
ual functional groups were found to concentrate primarily 
on regions with ETDs, surrounded by pristine (unfunc- 
tionalized) graphene areas. 

Given this current interest in defective graphene sys- 



tems, in the present work, we employ ab initio density 
functional theory (DFT) calculations to study the mor- 
phological evolution of the following graphene structures: 
(i) homogeneously-strained graphene sheets containing a 
Stone- Wales (SW) defectj^i and (ii) pristine and func- 
tionalized graphene sheets containing initial distributions 
of isolated topological defects (TD)j^® after introduction 
of holes and highly inhomogeneous bond deformations. 
We find that stress-relaxation of a graphene sheet con- 
taining SW defects may lead to the formation of ETDs, 
with morphological units that are very similar to those 
observed in the aforementioned experiments i^^— We un- 
cover the role of the SW defect as a stress-accumulation 
site that induces the bond rotation events which gen- 
erate the ETDs. Furthermore, we find that healing of 
small voids coupled with the relaxation of inhomoge- 
neously strained regions also lead to the formation of 
ETDs, revealing a rich variety of morphological patterns 
and plastic deformation mechanisms in pure and func- 
tionalized graphene structures. Our results indicate that 
the marked tendency of graphene to form TD clusters un- 
der strain is aided by functional groups, since the forma- 
tion energies of such defective structures are considerable 
smaller in graphane (a graphene sheet fully functional- 
ized with hydrogen atoms) and in graphenol (a graphene 
sheet fully functionalized with hydroxyl groups) than in 
graphene, in agreement with the experimental evidence 
in Ref. 5. 

Our calculations are performed using Kohn-Sham 
density functional theory;^ the generalized-gradient 
approximation's for the exchange-correlation functional, 
and norm-conserving Troullier-Martins pseudopoten- 
tials,= to describe the electron-ion interactions. We use 
the LCAO method implemented in the SIESTA code,— 
with a double-zeta pseudo-atomic basis set augmented 
with polarization orbitals, with an energy cutoff of 0.01 
Ry. Full structural optimization of atomic positions and 
supercell vectors is performed. For the relaxed struc- 
tures the total force on each atom is less than 0.02 eV/A 
and the pressure on the supercell is less than 1 kBar. 
The supercells employed are periodic along the graphene 



plane, being surrounded by a 30 A vacuum region in 
the transversal direction, with negligible interactions be- 
tween each layer and its periodic images. 



II. RESULTS AND DISCUSSION 

A. Homogeneous shear deformations of graphene 
containing Stone- Wales defects 

1. First principles results 

We consider initially a likely scenario, where homoge- 
neous shear deformations are imposed upon a graphene 
sheet containing a SW defect. SW defects are commonly 
found in graphene and nanotubes, and can be generated 
by irradiating the material with an 80 keV electronic 
beam, just below the knock-on threshold for sp^ bonded 
carbon materials.— i^>2i The impinging electrons give out 
the energy needed for the 90° rotation of a carbon-carbon 
(C-C) bond [as indicated in Figs.[lja) and (b)] that trans- 
forms four hexagons of the pristine graphene matrix into 
a SW defect. Since a SW is an incipient dislocation dipole 
in graphene (each pentagon-heptagon pair is one of two 
dislocations with opposite Burgers vectors comprising the 
SW defect), at first sight one expects shear deformations 
either to undo the C-C bond rotation and thus lead to 
self-annihilation of the two pentagon-heptagon pairs, or 
to lead to dissociation of the dislocation pair into its in- 
dividual components. 

In order to investigate this, we consider the two 



graphene supercells, each containing a SW defect, shown 
in Figllja) and (b). They differ by the orientation of the 
SW with respect to the x and y axes of the cell (indi- 
cated in the figure). We find that a 30% shear strain 
of the supercell shown in Fig. [ija) leads to the forma- 
tion of the ETD shown in Figs.[ll(c) and (d), containing 
clusters of fivefold and eightfold rings, while no morpho- 
logical transformation is observed when the same defor- 
mation is imposed on the supercell shown in Figlljb) . We 
shall point out that previous works have found pristine 
graphene to sustain deformations of up to 25% without 
yielding)^^"— a result that we also find in our calcula- 
tions, for graphene with SW defects. The morphological 
pattern observed in Fig. [D^c) is obtained by allowing the 
system to fully relax all the internal atomic coordinates, 
while maintaining the imposed homogeneous strain. No 
further changes in morphology are observed when full re- 
laxation of both atomic coordinates and strain (i.e. the 
supercell vectors) are allowed, as shown in Fig.[ljd). 

These results highlight a fundamental aspect of the 
plastic response of a graphene sheet to imposed stresses 
or strains: a tendency for the formation of extended topo- 
logical defects as a stress-relaxation mechanism, spawned 
by bond rotation events, in the above example with the 
formation of clusters of fivefold and eightfold rings. We 
note also that the transformation depends on the orien- 
tation of the shear deformation with respect to the pre- 
existing SW defect, that acts as stress- accumulation seed 
for the topological transformation. 



2. Analysis of elastic energy using a Keating model 




FIG. 1: Morphological transformation of a graphene sheet 
under strain, (a) and (b) show two images of graphene su- 
percells containing a Stone Wales defect. The bond rotation 
that generates the defect is indicated. The two cases differ 
by the orientation of the defect with respect to the x and y 
axes of the cell, (c) Transformation of the supercell in (a) 
with relaxation of internal atomic coordinates, under homo- 
geneous shear strain, (d) Fully-relaxed stress-free state of the 
structure in (c), after relaxation of cell vectors. 



The latter point can be illustrated by following the 
evolution of the elastic energy of the cell during the de- 
formation, as shown in Fig. [2] In the figure, we show a 
map of the elastic energy of the supercell at three dif- 
ferent stages of the plastic relaxation process, computed 
using a Keating model. ^^ Figure \^&) shows the elastic 
energy for the initial strained geometry. Note that, in 
this configuration, the strain energy is distributed nearly 
uniformly over the whole system, with a larger magni- 
tude at the SW defect. This pattern is preserved up 
to the point where a morphological transformation takes 
place. Figure EJb) shows the strain energy right after 
the transformation: the important feature here is that 
a change in the distribution of the elastic energy takes 
place, from nearly uniform before the transformation to 
a pattern where the energy is concentrated on the ensuing 
ETD, and where the bulk portion of the cell relaxes to a 
configuration of much smaller elastic strain. Relaxation 
of internal atomic coordinates, after the transformation, 
leads to further reduction of the strain in the cell. Fi- 
nally, in Fig. [5] (c) the imposed homogeneous strain is 
lifted and the system relaxes to a stress-free state, where 
the elastic energy in the bulk portions of the cell nearly 
vanishes, and most of the residual strain is located on the 
defect. 



B. Deformations of graphene containing small 
holes and topological defects 



1. Inhomogeneous deformations of graphene containing 
Stone- Wales defects 



We consider now more complex deformations of a 
graphene sheet, including small sized holes generated by 
partially ripping a graphene sheet, coupled with strongly 
inhomogeneous bond deformations. These can be consid- 
ered a model for severely strained and ripped graphene 
samples, such as those obtained by reduction of GO.^ 
In the following, we show that healing of such graphene 
fragments, combined with relaxation of the inhomoge- 
neous deformations, also lead to the formation of ETDs, 
in some cases with morphologies similar to those that 
have been observed in recent experiments. 



We start with the supercell shown in Fig. [IJb). By 
imposing a homogeneous shear strain, coupled with lines 
of strained C-C bonds (elongated by 27% to ^1.8 A) as 
indicated by dashed segments in Fig. [21Ja) and (c), and 
allowing for full relaxation of both internal coordinates 
and homogeneous strain, the system reaches a stress-free 
state where the dislocation dipole is dissociated into two 
individual dislocations separated by a lattice constant, as 
shown in Fig. ^h). When the strained bonds combine 
with a homogeneous shear oriented with respect to the 
SW defect as shown in FiglSlJc), the SW defect induces 
the nucleation of an ETD consisting of two side sharing 
octagons connected to two pairs of pentagons, each on 
one of two opposite sides of the octagon pair, shown in 
Fig.El^d). This ETD is similar in extent to those resulting 
from the 555-777 and 585 reconstructions of a divacancy 
in graphenei^ 




FIG. 2: Keating-energy analysis of stress relaxation by ex- 
tended topological defect formation. In the initial structure 
a homogeneous shear strain is imposed on a graphene cell 
containing a Stone- Wales defect, as shown in Fig. [T][a). The 
location of the defect is at the top, along the center of the 
cell, in the regions with higher elastic-energy histogram bars, 
(a) In the initial structure, the elastic-strain energy is nearly 
uniform in both defective and bulk regions of the cell, (b) 
Concentration of elastic energy, after morphological transfor- 
mation that generates an extended topological defect, in the 
process of stress relaxation under applied strain, (c) Relax- 
ation of bulk and defect strain in the stress-free metastable 
structure. Most of the residual strain is located on the defect. 



2. Inhomogeneous deformations of topological defect 
networks in graphene 

In order to further investigate the morphologies of de- 
formed graphene, and to access the role of functional- 
ization and preexisting topological defects, we choose 
as starting configurations the three topological-defect 
graphene networks^ shown on the left in Fig. 21 which 
were labelled Su, S31 and S22 in Ref. 18. We also 
consider the graphane-- (hydrogen-functionalized) and 
graphenol (hydroxyl-functionalized) versions of these 




FIG. 3: Combination of homogeneous and inhomogeneous de- 
formations of a graphene sheet: (a) and (c) show two images of 
homogeneously sheared supercells with a Stone Wales defect 
and inhomogeneous deformations generating lines of strained 
bonds, (b) Relaxation of the deformations in (a) leads to 
dislocation-dipole (Stone Wales defect) dissociation, (d) Re- 
laxation of the deformations in (c) leads to formation of an 
extended topological defect. 



structures, both at full coverage, with graphenol taken as 
a simple model for graphene oxide^^ A variety of defor- 
mation patterns, including small holes, as well as com- 
pressed and/or elongated bonds were initially imposed 
on these structures, and each system was allowed to relax 
(with optimization of both atomic coordinates and super- 
cell vectors) onto a stress-free metastable local minimum 
of the total-energy surface. 

Figure U shows starting configurations before and after 
introduction of holes and inhomogeneous compressions, 
and the resulting relaxed ETD configuration. In our 
calculations, these are, respectively, the lowest-energy 
ETD for each of the following three cases: (a) pristine 
graphene, (b) graphane, and (c) graphenol. The two 
functionalized cases in Figs.Sfb) and (c) show the forma- 
tion of a morphological pattern in the shape of a flower, 
where a square ring is surrounded by four heptagons, 
with pentagons connected at the perimeter. In the case 
of pristine graphene, the lowest-energy ETD we obtained 
in our calculations contains pentagon-heptagon clusters 
connected to pentagon-octagon units of the same mor- 
phology as the one that has been recently observed as 
the periodic unit of a domain-boundary in a graphene 
layer deposited on a nickel substratei^"— 

Further morphologies are obtained by considering the 
various deformation patterns described above. While the 
formation of the ETDs observed experimentally is largely 
dictated by kinetics, the thermodynamic stability of such 
structures can be analyzed by computing their forma- 
tion energies, Ef. In the following, we describe a to- 




FIG. 4: Lowest-energy extended topological defects obtained 
from regeneration of holes and response to inhomogeneous de- 
formations. Left column shows the parent geometry. Middle 
column shows imposed deformation. Right column shows re- 
sulting stress-free metastable ETD after full relaxation of in- 
ternal atomic coordinates and cell vectors, (a) Pristine (with- 
out functionalization) showing clusters of 5-7 pairs as well as 
a 8558 unit, (b) Hydrogen- functionalized and (c) Hydroxyl- 
functionalized graphene show a morphological pattern where 
a tetragon is surrounded by four heptagons, with pentagons 
connected to its perimeter. 



tal of 12 ETD configurations with Ef < 0.30 eV/atom. 
To put these formation energies in perspective, we quote 
the formation energy per atom of a fullerene molecule, 
0.40 eV/atom with respect to a graphene sheet, and 
the energy of the pentaheptite, a theoretical allotrope 
of graphene composed entirely of pentagons and hep- 
tagons, 0.24 eV/atom,^*''^'* both computed with the same 
methodology as the ETDs energies in this work. 



3. Formation energy and energetics 

In the present study, the formation energies per atom 
of the structures without functionalization are defined as 
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where Eetd is the ab initio total energy of the A^-atom 
supercell containing the ETD and ^c — Egraphene/N is 
the total energy per atom of an TV-atom bulk graphene 
cell. For systems with hydrogen functionalization, we 
take graphane as reference, such that the formation en- 
ergies per CH unit are defined as 
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where we keep fic as the energy per atom in graphene, 
and define fin = Kh - fJ-c, where ^icH = E graphane /N 
is the total energy per CH molecule for a bulk graphane 
calculation. For the hydroxyl-functionalized cases, we 
have: 
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for the energy per COH unit, where again we keep ^c 
as the energy per atom in graphene, and define ^oh = 
tJ-coH - Mc, where ^icOH = Egraphenoi/N is the total 
energy per COH radical in a bulk graphenol calculation. 

Table I shows the formation energies of the 12 ETDs, 
along with the energies of the parent geometries, S'12, 
531, and 5*22, both with and without functionalization. 
In some cases, regeneration of holes did not take place, 
a feature that is consistent with the presence of holes in 
the final state of the RGO samples in Ref. 5. These are 
indicated in Table I as unreconstructed (UNREC). 

Figure [5]Ja) shows the six structures derived from de- 
formed versions of the S12 geometry.— From the fig- 
ure, formation of pentagon-heptagons clusters is clearly 
the most common morphological patterns emerging 
from small-hole regeneration and relaxation of deformed 
graphene sheets, as found in recent experiments, but 
other morphologies are also manifest. We observe poly- 
gons ranging from tetragons to hendecagons, with color- 
coding indicated in the figure. The ETD in structure 
SA, with a nonagon connected to pentagons, is very sim- 
ilar to that present in the model of amorphous graphene 



recently proposed as a high-specific area material for su- 
percapacitor applications J5_ The five structures result- 
ing from deformations of the S31 geometry are shown in 
Fig-EJb). Linear polyacetylene-like chains are observed 
in the SG and SJ hydrogenated clusters, in structures 
with a strong sp'^-hibridization content, which are stabi- 
lized by the functional groups. Such polymer-like carbon 
chains were also obtained in molecular dynamics mod- 
els of RGO. Removal of functionalization in the SJ ge- 
ometry leads to the recovery of a threefold connected 
network of pentagon-heptagon clusters, in the form of a 
climbed dislocation dipole, of the type observed in irra- 
diated graphene samples^ Another striking geometry is 
the SI, which shows the presence of tetragons connected 
to decagons, with a reasonably low formation energy of 
0.3 eV/atom when functionalized with hydrogens. Note 
the presence of a network of large pores, which would 
make defective graphene sheets with such morphologies 
candidates for applications as supercapacitors^ and se- 
lective permeable membranesj^i^ 




FIG. 5; (a) Low-energy extended topological defects derived 
from deformed versions of the S12 structure (shown in the 
center). (b)Low-energy extended topological defects derived 
from deformed versions of the S31 structure (shown in the 
center at the top row). Color coding for topological defects, 
ranging from tetragons to hendecagons, is indicated. 



TABLE I: Formation energy Ef, in eV per formula unit, of 
extended topological defects in graphene. Eg is the elec- 
tronic band gap in eV. Structures without functionalization 
are metallic, except the S31. Functionalized structures are 
semiconductors or insulators, except the SE. 
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4-. Role of functionalization 



Regarding the role of functionalization, we observe 
in Table I that the formation energy decreases with 
functionalization in all cases, with the exception of the 
hydroxyl-functionalized S31, SD, and SI, which have 
slightly higher values of Ef than their pristine counter- 
parts. For the hydrogen-functionalized ETDs, an average 
reduction of 95 meV/formula-unit is observed in Table I, 
while for the hydroxyl-functionalized ETDs, we obtain 
an average reduction of 125 meV/formula-unit. In both 
cases rather sizable maximum reductions of Ef of 0.27 
and 0.37 cV/formula-unit occurs, respectively, for the 
H- and OH-functionalized SG morphologies, where func- 
tionalization leads to stabilization of polymer-like car- 
bon chains. Functionalization also changes the electronic 
structure: while all pristine ETDs are metallic, with the 
exception of the S31 parent geometry, the functionalized 
ETDs are all semiconductors or insulators, with the ex- 
ception of the SE. Table I shows the values of the gap for 
the semiconducting and insulating ETDs. This change 
in the electronic structure is connected with the incor- 
poration of an sp^ component in the electronic structure 
of these graphene materials, associated to the bonding of 
the functional groups to the graphene sheet. 



III. CONCLUSION 

In conclusion, ab initio calculations indicate that 
stress-relaxation of a graphene sheet containing topolog- 
ical defects may lead to the formation of extended topo- 
logical defects showing morphological units that have 
been observed in recent experiments. We uncover the 
role of topological defects as a stress-accumulation site 
that induces the bond rotation events that generate the 
ETDs. Furthermore, we find that healing of small voids 
coupled with the relaxation of inhomogeneously strained 
regions also lead to the formation of ETDs, revealing a 
rich variety of morphological patterns and plastic defor- 
mation mechanisms in pure and functionalized graphene 
structures. Our results indicate that the tendency of de- 



formed graphene sheets to form topological defect clus- 
ters is enhanced in the presence of functional groups, with 
a systematic reduction of the formation energies of such 
defective structures in graphane and in graphenol, when 
compared to graphene, in agreement with the available 
experimental evidence. 
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